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Light Gravitinos at Colliders and Implications for Cosmology 
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Light gravitinos, with mass in the eV to MeV range, are well-motivated in particle physics, but 
their status as dark-matter candidates is muddled by early-Universe uncertainties. We investigate 
how upcoming data from colliders may clarify this picture. Light gravitinos are produced primarily in 
the decays of the next-to- lightest supersymmetric particle, resulting in spectacular signals, including 
di-photons, delayed and non-pointing photons, kinked charged tracks, and heavy metastable charged 
particles. We find that the Tevatron with 20 fb~^ and the 7 TeV LHC with 1 fb~^ may both see 
evidence for hundreds of light-gravitino events. Remarkably, this collider data is also well suited to 
distinguish between currently viable light-gravitino scenarios, with striking implications for structure 
formation, inflation, and other early-Universe cosmology. 

PACS numbers: 14.80.Ly, 13.85.-t, 95.35.-|-d, 98.80.Cq 



I. INTRODUCTION 

Supersymmetry is one of the most promising ideas 
for new physics beyond the Standard Model. Super- 
symmetric theories that incorporate local supersym- 
metry (or supergravity) predict the existence of the 
gravitino, the spin-3/2 superpartner of the graviton. 
When supersymmetry is broken, the gravitino acquires 
a mass through the super-Higgs mechanism, "eating" 
the spin-1/2 goldstino, the Goldstone fermion associ- 
ated with spontaneously-broken local supersymmetry [J- 
0]. In contrast to other superpartners, the gravitino 
can have a mass m^, that is not at the weak scale 
Wwoak ~ 100 GeV — 1 TeV, and viable models exist for 
gravitino masses as low as the eV scale and as high as 
100 TeV. In this work, we consider light gravitinos, with 
mass in the eV to MeV range. Such gravitinos are highly 
motivated in particle physics, as they emerge in models 
with gauge- mediated supersymmetry breaking (GMSB), 
in which constraints on flavor violation are naturally sat- 
isfied [i-[i3i. 

Light gravitinos also have cosmological motivations. In 
particular, they are the original supersymmetric dark- 
matter candidate [Tlj . Assuming a high reheating tem- 
perature, gravitinos are initially in thermal equilibrium 
and then freeze out while still relativistic. As we discuss 
in detail below, their resulting relic density is 
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where g*s,f is the number of relativistic degrees of free- 
dom at freeze out, and has been normalized to the total 
number of degrees of freedom in the Standard Model. 
When originally proposed in the 1980's, uncertainties in 
h and the total matter relic density allowed ^ keV. 
This led to a simple and attractive gravitino-dark-matter 
scenario, consistent with standard Big Bang cosmology, 
in which the Universe cooled from some high tempera- 
ture, and keV gravitinos froze out and now form all of 
the dark matter. 



In the intervening years, however, a variety of astro- 
physical constraints have greatly complicated this pic- 
ture. First, the dark-matter relic density is now known 
to be iluuh'^ — 0.11. Second, constraints on structure 
formation, as probed by galaxy surveys and Lyman-a 
forest observations, require that the bulk of dark matter 
be cold or warm [l^l . As we will discuss more fully below, 
this leads to three scenarios of interest: 

1- "^G ^ 15 — 30 eV: Gravitinos are produced by 
the standard cosmology leading to Eq. ([T]); they 
are hot dark matter, but their contribution is small 
enough to be consistent with the observed small- 
scale structure. Some other dark-matter particle is 
required. 

2. 15 — 30 eV < rriQ < few keV: Non-standard cos- 
mology and a non-standard gravitino production 
mechanism are required, both to avoid overclosure 
and to cool the gravitinos to satisfy small-scale- 
structure constraints. Some other dark-matter par- 
ticle may be required. 

3- "T-Q > few keV: Non-standard cosmology is re- 
quired to dilute the thermal relic density of Eq. ([T]) . 
Gravitinos produced by thermal freeze out are cold 
enough to be all of the dark matter. 

Note that the original "keV gravitino" scenario, previ- 
ously favored, is now the most disfavored, in the sense 
that it is excluded by both overclosure and small-scale- 
structure constraints. All of the possibilities are rather 
complicated, however, as in each case, some additional 
physics is required, either to provide the rest of the dark 
matter or to modify the history of the early Universe to 
allow gravitinos to be all of the dark matter. 

In this paper, we discuss how collider data may help 
clarify this picture. Light gravitinos are primarily pro- 
duced at colliders in the decays of the next-to-lightest 
supersymmetric particle (NLSP). It is a remarkable co- 
incidence that modern particle detectors, with compo- 
nents placed between 1 cm to 10 m from the beamline, 
are beautifully suited to distinguish between the NLSP 
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decay lengths predicted in scenarios 1,2, and 3. For ex- 
ample, the decay length of a Bino NLSP decaying to a 
gravitino is [13] 



scale of the superpartner masses, is [22 



CT ~ 23 cm 



"G 



100 cV 



100 GeV 



(2) 



This implies that scenarios 1,2, and 3 make distinct pre- 
dictions for collider phenomenology, and the identifica- 
tion of the gravitino collider signatures realized in nature 
may have far-reaching implications for the early Universe. 

Of course, this requires that gravitinos can be pro- 
duced in sufhcient numbers and distinguished from Stan- 
dard Model backgrounds. In this work, we determine 
event rates for a variety of signatures, including prompt 
di-photons and delayed and non-pointing photons (rel- 
evant for neutralino-NLSP scenarios), as well as kinked 
charged tracks and heavy metastable charged particles 
(relevant for stau-NLSP scenarios). We present results 
for an assumed final Tevatron dataset (20 fb~^ of 2 TeV 
pp collisions), an early LHC dataset (1 fb"""^ of 7 TeV pp 
collisions), and a future LHC dataset (10 fb^^ of 14 TeV 
pp collisions). We find that the final Tevatron and early 
LHC data have roughly equivalent sensitivity to these 
events, with both capable of seeing hundreds of distinc- 
tive light gravitino events. The full LHC data greatly ex- 
tends the reach in parameter space, and may also allow 
precision measurements of NLSP lifetimes and gravitino 
masses. 

We begin in Sec. |TT] by reviewing the cosmological 
bounds on light gravitinos and discussing how these 
bounds are relaxed in early-Universe scenarios that dif- 
fer from the canonical one. In Sec. IIIII we then discuss 
NLSP decays to gravitinos, GMSB models, and current 
collider constraints. In Sec.|TV]we present our results for 
the number of light-gravitino events at colliders, based 
on collider simulations, and discuss the cosmological im- 
plications. We summarize our conclusions in Sec. [V] 



II. LIGHT GRAVITINO COSMOLOGY 
A. Canonical Scenario 

1. Relic Abundance 

In the currently canonical scenario, after inflation, the 
Universe is reheated to a temperature that is as- 
sumed to be far higher (e.g., 10^^ or 10^^ GeV) than 
the weak scale. During this phase, inelastic scattering 
processes and decays can convert Standard-Model par- 
ticles in the thermal bath into gravitinos [iHl3. The 
rate Cq per unit volume for production of light graviti- 
nos (strictly speaking, only the spin-1/2 goldstino compo- 
nents) can be calculated by considering all such processes, 
which primarily involve strong [l^ and electroweak gauge 
bosons 20j, |2l| , as well as top quarks [12] . The total re- 
sult, valid in the limit T ^ ?tt-susy, where msuSY is the 
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where Mpi ~ 1.2 x 10^^ GeV is the Planck mass. Here 
we have assumed that the gaugino masses Wg,!, '7ig^2, 
and TOg^3 and the tri-linear scalar coupling At are at a 
common mass scale. For simplicity, we have set them 
equal to a universal gaugino mass rrig. 

The evolution of the gravitino number density via 
these production processes, and their inverses, is gov- 
erned by the Boltzmann equation 



driQ 
~~dt 



ZHtlq — Cq — Fn^ 



(4) 



where H is the Hubble expansion rate and F is the rate 
of processes that annihilate gravitinos. The 3Hn^ term 
accounts for dilution of the number density due to cosmo- 
logical expansion. If F 3> H , gravitinos are in thermal 



equilibrium, Tn^ 



and their number density (the 



solution to the Boltzmann equation) is 
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^^r3~0.24T3 
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Here we used 5 = 2, since it is primarily the spin-1/2 
goldstino components that are produced thermally. 

The rate F at which a given gravitino is destroyed in 
the plasma is then 



(6) 



Since F a and F a T^, the ratio T/H (xT is largest 
at the highest temperatures. Thus, if F(Tr) > H{Tr) at 
reheating, then gravitinos come into thermal equilibrium 
shortly after reheating. During this era, the expansion 
rate is given by i? ~ 1.66 ^y^T^/Mpi; assuming reheat- 
ing temperatures Tr ^ TeV, at which all particles in the 
minimal supersymmetric Standard Model (MSSM) are 
relativistic, we set the number g^, of relativistic degrees 
of freedom to .g*(rfl) ~ 228.75. Comparing T{Tr) and 
H{Tb), we then see that if the reheating temperature 
satisfies 



Tr > Tf = 5 GeV 



keV 



TeV 



n 2 



(7) 



then gravitinos come into thermal equilibrium after re- 
heating. Recalling that the production rate, Eq. (jSj, used 
here is valid only for T ^ tonlsp (i.e., T ^ 10 TeV), we 
conclude that for weak-scale gluino masses, light graviti- 
nos with vtiq < MeV will come into thermal equilibrium 
if the reheating temperature is Tr ^ 5 x 10^ GeV. 

The creation/annihilation rates for gravitinos at tem- 
peratures T < 10 TeV have not yet been calculated, 
and so the precise temperature at which gravitinos 
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freeze out (which occurs when F ~ H) cannot yet be 
determined. Still, gravitinos are produced and de- 
stroyed individually, requiring (from R-parity conserva- 
tion) that each creation/destruction is accompanied by 
creation/destruction of some other supersymmetric par- 
ticle. Therefore, the freezeout temperature Tf cannot be 
much lower than the mass jtinlsp of the NLSP, as the 
equilibrium abundance of SUSY particles then decreases 
exponentially. We thus conclude that the freezeout tem- 
perature for light gravitinos falls roughly in the range 
10 GeV <Tf < 10 TeV. 

With this range of freezeout temperatures, Tf ^ m^, 
so gravitinos are relativistic when they freeze out. The 
relic gravitino density is then [ll| 



0.1 



100 cV 



106.75 

9*s,f 
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the standard result for hot relics with g = 2, where g*sj 
is the number of relativistic degrees of freedom when 
the gravitinos freeze out. If gravitinos freeze out when 
all of the MSSM degrees of freedom are relativistic, 
Tf ^ TOsuSY, and g*s,f = 228.75. However, it is more 
likely that freeze out occurs at T ~ wnlsp ^ 100 GeV, 
when g^sj ~ 100 [H. 



2. Cosmological Constraints 

Given that current cosmic-microwave-background 
(CMB) and structure-formation measurements constrain 
the density of dark matter to be JIdm^^ — 0.11, Eq. 
implies an upper bound rriQ < 200 eV. The upper 
limit uiq ~ 200 eV is saturated if the gravitino makes 
up all of the dark matter and freezeout occurs when 
9*sj = 228.75 is the maximum value allowed in the 
MSSM. However, a gravitino of this mass would be hot 
dark matter. It would smooth density perturbations on 
scales probed by galaxy surveys and the Lyman-a for- 
est to a degree that is highly inconsistent with data. As 
is reduced from this upper limit, the smoothing scale 
is increased (the gravitinos get "hotter"), but the grav- 
itino abundance is reduced, thus making the magnitude 
of the smoothing smaller. A combination of data from 
the cosmic microwave background, galaxy surveys, and 
the Lyman-a forest constrain the contribution of a hot 
component of dark matter to be < 15% [l^l, implying 
for g*s,f — 100 that ttIq < 15 eV. This suggests that the 
most conservative upper bound is given by ttiq < 30 eV, 
in the case that g*s.f — 200 approaches the maximal 
value allowed in the MSSM. Therefore, in this canonical 
scenario, thermal gravitinos with mass less than 30 eV 
make up only a fraction of the dark matter, thus requir- 
ing some other particle to be the cold dark matter. This 
is the first cosmological scenario listed in Sec. HI 

We conclude by noting that future astrophysical data 
are likely to improve. And while the current sensitivity 
is to gravitino masses as small as 15 — 30 eV, it is fore- 
cast that next-generation experiments may be sensitive 



to gravitino masses as small as 1 eV [2J|. A collider de- 
tection of a gravitino in the mass range m^^ = 1 — 30 eV 
would thus lead to testable consequences in forthcoming 
cosmological data. 



B. Non-standard Early-Universe Scenarios 

There are several ways in which the early- Universe pro- 
duction of relic gravitinos could differ from the canonical 
scenario outlined above. Thus, there are scenarios in 
which a gravitino of mass ttiq > 30 eV, ruled out in the 
canonical model, could be cosmologically consistent or, 
better yet, completely compose the dark matter. 

Let us first consider scenarios in which the gravitinos 
reach thermal equilibrium in the early Universe, since 
most observational constraints are strictly valid only un- 
der this assumption. As mentioned above, if we only 
consider particles in the MSSM, then g*sj ^ 228.75, and 
Eq. ^ suggests an upper limit of itiq < 200 eV from the 
relic abundance constraint. One way to evade this limit 
is to simply consider higher values of g*sjj i-e., graviti- 
nos decouple and freeze out earlier than in the canonical 
scenario. This may be possible in models with more de- 
grees of freedom than the MSSM. More massive graviti- 
nos that decouple earlier may then be viable, if they have 
an abundance that obeys the constraint JIdm^^ ^ 0.11. 

It is possible that this constraint is saturated and that 
these heavier gravitinos entirely compose the dark mat- 
ter. Of course, we must still require that these heavier 
gravitinos are not so hot as to erase structure to a de- 
gree that contradicts observations. The same combina- 
tion of CMB, galaxy-survey, and Lyman-a-forest data 
that was used to constrain itIq < 30 eV in the canonical 
scenario can also be used to constrain the gravitino mass 
in this early-decoupling scenario, assuming that thermal 
gravitinos make up all of the dark matter. With this 



assumption, Ref. 



find 



> 550 eV, using a selec- 



tion of Lyman-a data. The same authors later find a 
stronger constraint, > 2 keV, with SDSS Lyman-a 
data |25l . [26| , a result slightly weaker than a bound on 
warm-dark-matter models obtained by Ref. [27| . A num- 
ber of other small-scale observations also seem to support 
that rriQ > few keV under these assumptions [28| . 

We thus conclude that if 30 eV < ttiq < few keV, then 
thermal gravitinos are too warm to be the only com- 
ponent of the dark matter, regardless of whether or not 
they have the correct abundance. Gravitinos in this mass 
range would only be viable if some other non-standard 
early-Universe process cools them, or if there is an addi- 
tional cold component. This is the second scenario men- 
tioned in SecHl However, if rriQ > few keV, then graviti- 
nos may be sufficiently cold, and may in early-decoupling 
scenarios have the right abundance, to be the dark mat- 
ter. This is the third scenario outlined in Sec. HI 

Of course, aside from early decoupling, there are other 
non-standard mechanisms that can reduce the gravitino 
abundance. For example, recall that we have no empiri- 
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cal constraints to the early Universe prior to the epoch of 
big-bang nucleosynthesis (BBN), at which T > few MeV 
[29] ■ Thus, some entropy-producing process prior to BBN 
could also dilute the gravitino abundance. It is possi- 
ble that there may be some exotic early- Universe physics 
that conspires to produce the same effect. One relatively 
simple possibility is that the reheating temperature is 
low. If the reheating temperature is smaller than the 
freezeout temperature, then gravitinos will never come 
into thermal equilibrium, and their relic abundance will 
thus be accordingly smaller [30l - l36j . The only catch is 
that for the light gravitinos we consider here, the reheat- 
ing temperature must be unusually low for this to oc- 
cur. For example, if rriQ = keV and rug = 300 GeV, then 
Eq. ([T]) suggests that the reheating temperature must be 
Tji < 50 GeV. However, recall that this estimate may 
not be strictly valid at T < 10 TeV, as we have already 
noted. Thus, a more careful calculation of the production 
rate of light gravitinos at low reheating temperatures may 
be necessary. Nevertheless, such low reheating tempera- 
tures have been considered [s^, and Ref. [11] has exam- 
ined an explicit low-reheat scenario in which a gravitino 
of mass ttiq = 1 — 15 keV can have the right abundance 
to be the dark matter. 

Finally, we also note that there may be additional 
mechanisms affecting the generation of gravitinos. For 
example, in our discussion we have neglected the non- 
thermal contribution to the gravitino abundance from 
out-of-equilibrium decays of other supersymmetric parti- 
cles. There may also be other significant modes of grav- 
itino production or dilution, including processes involv- 
ing the messenger particles responsible for GMSB [sol - 
|44| . nonthermal production via oscillations of the infla- 
ton field [i^, and various other mechanisms [4^. There 
may thus be other reasons why the gravitino abun- 
dance or temperature differs from those in the canoni- 
cal thermal-production scenario; this may be true even if 
TOg < 30 eV. 

To summarize, in the canonical model, gravitinos are 
required to have mass < 30 eV and form only a frac- 
tion of the dark matter. Gravitinos with mass range 
TOq > 30 eV would require non-standard physics or cos- 
mology to reduce their abundance or temperature to 
agree with observations. Below we discuss collider sig- 
natures of light gravitinos. We close here by noting that 
such collider data may, if gravitinos are discovered, thus 
help discriminate between the diversity of early-Universe 
scenarios for gravitino production. 



III. LIGHT GRAVITINOS AT COLLIDERS 



sion 2, as 



\/3M, 



240 eV 



103 TeV 



(9) 



Stt ~ 2.4 X 10^* GeV is the reduced 



where M* = Mpi/\ 
Planck mass. 

The interactions of weak-scale gravitinos are of gravi- 
tational strength, as expected since they are the super- 
partners of gravitons. However, the couplings of the gold- 
stino are proportional to 1/i^ [13, \^ ■ The interactions 
of light gravitinos are therefore dominated by their gold- 
stino components, and may be much stronger than grav- 
itational. Decays to gravitinos are faster for light grav- 
itinos. 

For reasons to be discussed below, we will focus on 
cases where the NLSP is either the neutralino or the stau. 
For a neutralino NLSP that is dominantly a Bino, the 
decay widths to gravitinos are [H, 

cos^ 9v 



T{B 7G) 



r{B ZG) = 



167rF2 
sin^ 9wm^l 
167rF2 



1 - 



(10) 
(11) 



where 9w is the weak mixing angle. For rrij^ < mz, de- 
cays to Z bosons are negligible or kinematically forbid- 
den, and the corresponding decay length is 



cr ~ 23 cm 



r "^G 1 


2 


'100 GeV" 


LlOO evJ 







(12) 



For heavier neutralinos, the Z mode may be signifi- 
cant; for very heavy Binos, the branching ratio for this 
mode is B{Z) ~ sin^ 9w — 0.23. Decays to l~^l~G, where 
I = (e,/i,orT) is a charged lepton, and hG may also be 
possible; however, these modes have branching ratios of 
^ 0.01 and ~ 10~^, respectively. 

For stau NLSPs, the decay width is |50i] 



T{t tG) 



corresponding to a decay length 



cr ~ 18 cm 



r '"g 1 


2 


'100 GeV" 


Lioo cvJ 
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(13) 



(14) 



As anticipated, in both the neutralino- NLSP and stau- 
NLSP scenarios, the decay lengths for gravitinos in the 
cosmologically interesting range correspond to distances 
that bracket the size of collider detectors. 



Mass and Interactions 



B. GMSB Models 



The gravitino mass is determined by the super-Higgs 
mechanism. In simple models, it is given in terms of the 
supersymmetry-breaking scale F, which has mass dimen- 



Light gravitinos are expected to be dominantly pro- 
duced at colliders in the cascade decays of strongly- 
interacting superpartners, such as squarks and gluinos. 
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Collider constraints therefore depend on the full super- 
partner spectrum, and so are mo del- dependent. Follow- 
ing most of the literature, we will work in the framework 
of minimal GMSB, and so we briefly review its features 
here. 

Typical GMSB models are characterized by a hid- 
den sector, a messenger sector, and a visible sector, 
the MSSM. Supersymmetry breaking is triggered by a 
hidden-sector gauge-singlet superfield S acquiring the 
vacuum expectation value S = M + O^Fg. This then 
generates masses for the messenger-sector fields M„iess = 
AM, where A is a coupling in the superpotential. These 
in turn generate masses for the visible-sector superpart- 
ners that are roughly a loop factor times A = Fs/M, and 
so A 100 TeV. Note that Mmcss > A is generally as- 
sumed. 

In the minimal GMSB framework, the entire super- 
partner spectrum is specified by the parameters 

A, Mmcss, N5, tan/3, sgn(/x), Cg^av ■ (15) 

Here, A and Mmcss are as described above; masses and 
couplings are generated at Mmcss and then evolved to the 
weak scale via the renormalization group. The number of 
messenger superfields is given by N5 , the effective number 
of 5 -I- 5 representations of SU(5). The Higgs sector is 
specified by the usual parameters tan/3 and sgn(/i). The 
last parameter is 

^ (16) 



-grav 



where F — {Fg + F^)^/"^ is the total supersymmetry- 
breaking vacuum expectation value, which appears in 
Eq. ([9]). These relations imply 

(17) 



We expect Cg^av ^ 1, since F > Fs and A < 1, and in the 
minimal case that there is only one non-zero i^-term, we 
expect Cgiav ~ 1- 

The superpartner masses are determined by the pa- 
rameters of Eq. (fT5)) : for details, see Ref. [5l|. Here we 
note only two things. First, the superpartner masses are 
determined by g aug e couplings. Thus, although, for ex- 
ample, chargino [52| | and sneutrino [s^ NLSPs have been 
considered, the canonical NLSP candidates are those 
with only hypercharge interactions, namely, the Bino and 
right-handed sleptons. Among the right-handed slep- 
tons, the stau is typically the lightest, as renormalization- 
group evolution and left- right-mixing effects both de- 
crease the stau mass relative to the selectron and smuon, 
and so we will focus on the Bino-NLSP and stau-NLSP 
scenarios Q Second, the Bino and stau masses are pro- 



^ Note that in the "slepton co-NLSP" scenario, where the three 
charged sleptons are degenerate to within the mass of the tau, 
the number of e — > eG and /j — > decays may be comparable 
to that of the f — )■ tG decay that usually dominates gravitino 
production. 



portional to iVs and respectively. For iVs = 1, the 

NLSP is the Bino in minimal GMSB, but for iVs > 1, the 
stau may also be the NLSP; see, for example. Fig. 1 of 
Ref. [53. 

Thus, to study the Bino- like neutralino-NLSP scenario, 
we will choose — 1] likewise, we choose — A to 
study the stau-NLSP scenario. For both scenarios, we 



fix tan/3 = 20, > 0, and Cg 



1. We let A and M„ 



be free parameters. Note that the overall mass scale of 
the supersymmetric partners is roughly proportional to 
A, while the gravitino mass depends on both A and Afmcss 
as in Eq. ([TT]). Thus, scanning over the free GMSB pa- 
rameters will allow us to explore collider signals for a 
range of masses. We shall now examine the existing col- 
lider constraints on the parameter spaces of these two 
scenarios. 



C. Current Collider Constraints 

The high-energy collider signals of GMSB and graviti- 
nos are well studied |49l|50l . [5J-[63| ; for a review of cur- 
rent bounds, see Ref. [68|. Here we summarize the most 
relevant results for the models and signals we consider 
below. 

We shall discuss GMSB signals in more detail below, 
but we summarize them briefly here. In the neutralino- 
NLSP scenario, there are several possible signals. For 
short-lived neutralinos, all supersymmetry events include 
two prompt high energy photons. For longer-lived neu- 
tralinos that travel a macroscopic distance before decay- 
ing to photons in the detector, delayed or non-prompt 
photons are possible. The stau-NLSP scenario may also 
lead to a variety of signatures, depending on the stau 
lifetime, including acoplanar leptons, tracks with large 
impact parameters, kinked charged tracks, and heavy 
metastable charged particles. 

Several studies have attempted to place constraints 
on GMSB models by searching for these signals. Given 
that we will scan over a large range of the GMSB pa- 
rameter space, we are primarily interested in constraints 
that are generally valid over this entire range. We shall 
thus focus on limits from LEP studies, based on an in- 
tegrated luminosity of 628 pb~^ at center-of-mass ener- 
gies of 189 — 209 GeV, which combined searches for both 
GMSB and neutral- Higgs signals [HllO]. The relevant 
results for our models are the lower limits of A > 70 TeV 
for our neutralino-NLSP model, and A ^20 TeV for our 
stau-NLSP model; see Fig. 6 of Ref. [60|. These con- 
straints on A are valid for all values of Mmess we include 
in our scan. Therefore, the allowed region of Mmess — A 
parameter space is constrained by these LEP bounds. 

However, there are also a number of studies that fo- 
cused on constraining specific benchmark models [Tlj . 
which occupy certain points or lines in the GMSB pa- 
rameter space. Although these constraints cannot be di- 
rectly applied to our models, we discuss them to get an 
idea of the robustness of the LEP bounds on our param- 
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eter space. 

Of these benchmark-model constraints, the best col- 
Uder bounds on di-photon events are from the Tevatron, 
inchiding a DO search based on an integrated luminosity 
of 1.1 fb"^ [72] and a CDF search based on 2.6 fb"^ (tJ- 
The DO and CDF bounds, when interpreted assuming the 
benchmark GMSB model SPS 8 from Ref. [ll|, lead to 



lower bounds on the Bino mass of 125 GeV and 150 GeV, 
respectively. For longer-lived neutralinos that travel a 
macroscopic distance before decaying to photons in the 
detector, a CDF search for delayed photons, based on 
570 pb~^ of data, established lower bounds on from 
70 to 100 GeV for neutralino decay lengths between 20 
cm and 6 m, again when interpreted in the context of 
SPS 8 [zl- 

Searches for heavy metastable charged particles have 
also been performed at DO [t^ , assuming the benchmark 



GMSB model SPS 7 from Ref. |7l|. A similar search 
was performed at CDF 76] , but did not interpret results 
in the context of GMSB models. Based on ^ 1 fb^^ of 
data, and assuming only Drell-Yan slepton production, 
the constraints resulting from these two searches are not 
competitive with the LEP bounds stated previously. 

Thus, we shall take the more general LEP bounds as 
constraints on the two models we consider in this work, 
and shall further take only conservative values of the 
lower limits. For the neutralino NLSP model, we shall 
only scan the parameter space with A > 80 TeV, which 
should be comfortably allowed by the LEP bounds. How- 
ever, we acknowledge that it is possible that the Tevatron 
data may exclude a small range of NLSP masses within 
this parameter space comparable to that ruled out in the 
benchmark model (i.e., < 150 GeV), should this data 
be reanalyzed in the context of our models. For the stau 
NLSP model, we shall scan over A > 30 TeV. Given that 
the current Tevatron constraints are not competitive with 
the LEP bounds, all of this parameter space should be 
allowed. As we will see, hadron colliders have bright 
prospects for probing the parameter spaces of these mod- 
els. 



IV. TEVATRON AND LHC PROSPECTS 

A. Gravitino Signals 

The collider signal of a supersymmetric particle de- 
caying to a gravitino can be classified by (1) the distance 
from the interaction point at which the decay occurs, and 
(2) the nature of the accompanying Standard Model de- 
cay products. The former is determined by the gravitino 
mass and the masses of the decaying supersymmetric par- 
ticles, as well as the speed with which the decaying parti- 
cles are produced. The latter is determined primarily by 
the nature of the NLSP. We shall define and investigate 
the following categories of events: 

1. Prompt di-photons (in neutralino NLSP models): 
Events in which two photons are produced (via a 



pair of neutralino decays to gravitinos) within dp^ 
of the interaction point. We take dpi- — 1 cm as a 
conservative estimate of the distance to which the 
origin of any photon can be resolved in detectors 
at the Tevatron and the LHC. Note that here and 
below, we cut on the total distance traveled by the 
NLSP before it decays, not its (transverse) distance 
from the beamline when it decays. 

2. Non-prompt photons (neutralino NLSP): Events in 
which at least one photon is produced at a mid- 
detector distance ddccay away from the interaction 
point, where dpr < ddccay < dnp, and dnp is the 
maximum distance from the interaction point at 
which a photon can be observed. We conserva- 
tively take dnp = 3 m, roughly the outer radius 
of the hadronic calorimeters at both the Tevatron 
and the LHC. (Note that although the calorimeters 
in the ATLAS detector at the LHC actually extend 
to ~ 4 m, those in the CMS detector only extend to 
~ 3 m; we have thus taken the more conservative 
3 m as our cut.) Photons may also convert and be 
seen in the muon chambers, extending the sensitiv- 
ity to decays '-^ 10 m from the interaction point, 
but we neglect this possibility here. Here we also 
take dpr = 1 cm. 

Note that this category of events encompasses 
both non-pointing photons and delayed photons. A 
non-pointing photon is simply a photon that does 
not spatially point back to the interaction point. 
A delayed photon has the further distinction of be- 
ing produced only after a significant temporal delay 
following the time of the initial collision. This may 
occur when the particle that decays to the photon 
is produced with a low speed, so that it takes a non- 
negligible amount of time to travel away from the 
interaction point before it decays. If this amount 
of time is comparable to the time between collision 
events, it may be difficult to properly identify the 
delayed photon with its originating event. 

3. Non-prompt leptons (stau NLSP): Events in which 
at least one charged lepton is produced (via 
charged-slepton decays to gravitinos) at a mid- 
detector distance ddccay, where dpi- < ddccay < ^np 
as before. We take dnp = 5 m and d^p = 7 m as the 
outer radii of the muon chambers in the detectors 
at the Tevatron and the LHC, respectively. (As 
above, although the muon chambers in the ATLAS 
detector at the LHC extend to ^ 10 m, those in 
the CMS detector only extend to ~ 7 m; we take 
the more conservative 7 m as our cut.) We again 
take dpr = 1 cm. Each of these events produces a 
distinctive charged track with a kink due to the 
momentum carried away by the gravitino. 

As above, both non-pointing and delayed 
events are included in this category. Furthermore, 
we include all generations (e, fi, and r). As men- 
tioned previously, the stau is generally the lightest 
slepton, and hence we expect the majority of the 
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decays in the stau-NLSP scenario to be of the form 
f — >■ tG. Ahhough the heavier sleptons I — {eor p) 
may also decay to ZG, the branching ratio of this 
decay is generally suppressed compared to the de- 
cay to a lepton and a neutralino, i.e., ly^ . If the 
latter is kinematically forbidden, then the 3-body 
decays to It~t^ or Zr+f~ dominate instead. How- 
ever, as the mass splitting between the stau NLSP 
and the heavier sleptons decreases, these 3-body 
decays become less dominant (becoming kinemat- 
ically forbidden if the mass splitting becomes less 
than the tau mass). The decays to IG may then 
occur if the heavier sleptons I are produced at the 
end of a decay chain. 
4. Metastable sleptons (stau NLSP): Events in which 
at least one charged slepton passes through the en- 
tire detector before decaying to a charged lepton 
and a gravitino. That is, the gravitino is produced 
at ddecay > c?ms, whcrc dms IS the distance to the 
outer edge of the detector. We take dms = 5 m and 
dms = 10 m as conservative estimates of the sizes of 
the detectors at the Tevatron and the LHC, respec- 
tively. All generations (e, jl, and f ) are included. 
These events will produce charged tracks with a 
relatively large radius of curvature. 

For this category, we impose a further cut, re- 
quiring that the speeds /3 of the sleptons satisfy 
the criteria /3iower < < /?uppcr- The lower cut re- 
moves slower sleptons, which may be identified with 
the incorrect collision event. The higher cut re- 
moves faster sleptons, which may be misidentified 
as muons. We take typical values /3iowcr = 0.6 and 
/Supper = 0.8. Note, however, that Ref. |67| suggests 
a new search strategy that may be sensitive to even 
higher values of (3. 



All of these events will also be distinguished by miss- 
ing energy and momentum carried away by the graviti- 
nos. Note that these categories are chosen to be illus- 
trative of the variety of signals that may be observed, 
and that they are not comprehensive — we do not inves- 
tigate prompt di-lepton events or neutralino decays to Z 
bosons, for example. Furthermore, the categories are not 
mutually exclusive; for example, one may easily have a 
single event in which both a non-prompt lepton and a 
metastable slepton are produced. It is also clear that the 
relevant detector systematics and backgrounds will also 
be different for each category. 

This categorization of events is somewhat oversimpli- 
fied, as it is based primarily on cuts on the decay length. 
Certainly, additional cuts will be required in a realistic 
analysis, possibly reducing the number of detected sig- 
nals. However, we shall soon see that these simple cat- 
egories align with the three cosmological scenarios out- 
lined previously. 



B. GMSB Scan and Collider Simulations 

We now calculate the event rates for these gravitino 
signals in a parametrized GMSB model. A large number 
of programs have been written for the numerical com- 
putation of the mass spectra and collider predictions for 
parametrized supersymmetric models l77l 78| . In this 
paper, we use ISAJET 7.80/ISASUSY [zg to generate 
mass spectra and decay branching ratio tables. ISASUSY 
properly includes a number of 3-body decay processes 
relevant for gravitino phenomenology that are missing in 
other branching ratio programs. 

IS A JET/IS ASUS Y takes values of the GMSB param- 
eters listed in Eq. (llSp as input. As discussed previously, 
here we focus on parameterizations that fix a subset of 
the GMSB parameters, resulting in either a neutralino 
or a stau NLSP. We then scan over Mmess and A (requir- 
ing that Mmess > A), resulting in spectra with a range 
of gravitino and NLSP masses. The correspondence be- 
tween the Mmess — A scan and the resulting — tonlsp 
parameter space is shown in Fig. [TJ 

We then take the spectra and decay tables output by 
ISAJET/ISASUSY and use them as input for the Monte 
Carlo event generator PYTHIA 6.4.22 [8C|, including all 
supersymmetric processes available therein. For a given 
center-of-mass energy, PYTHIA can simulate a given 
number of collision events, giving a complete record of 
the various decay chains and final products generated in 
each event and an estimation of the various production 
cross sections. From this record, we can identify the su- 
persymmetric "mother" particles that decay to directly 
produce gravitino and Standard Model "daughter" par- 
ticles in each individual event. We can also find the de- 
cay length ddccay away from the interaction point that 
each mother particle travels before decaying to produce 
a gravitino. Thus, for any number of simulated events, 
we can find the fraction that fall into each of the above 
categories. The expected number of signals from each 
category is then given by the respective fraction multi- 
plied by the total number of supersymmetric events. We 
can also calculate the average (ddccay) of the decay length, 
taken over all supersymmetric events. 



C. Cosmological Implications 

The results of the scan are shown in Figs. [2][5] We 
can see that the simple categorization of collider sig- 
nals by decay-length cuts corresponds surprisingly well 
with the categorization of cosmological scenarios out- 
lined previously. For example. Fig. [2] shows that the 
observation of hundreds of prompt events suggests that 
the first cosmological scenario {m^ < 30 eV) is likely 
to be valid. Likewise, the second cosmological sce- 
nario (30 eV < rriQ < few keV) will be implied by the 
observation of a large number of non-prompt events, as 
demonstrated by Figs. [3] and [H Finally, that the ob- 
servation of a large number of metastable sleptons sup- 
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Neutralino NLSP Stau NLSP 




mg (keV) mg (keV) 

FIG. 1. Plots showing the mapping between the — rriNLSP and the Mmoss — A GMSB parameter spaces, for the neutralino- 
NLSP scenario with — 1 (left) and stau-NLSP scenario with A^s = 4 (right), where we fix tan /? — 20, A* > 0, and Cgrav = 1 in 
both cases. Contours of constant Mmess/A (dashed black) and A (dash-dotted black) are shown. The region in the upper-left 
corner is disallowed by theory, while the region at the bottom is excluded by experiment (using the conservative constraints 
mentioned in the text). 



ports the third cosmological scenario (m^j > few keV) 
can be seen in Fig. [S] We emphasize that this correspon- 
dence is not strongly dependent on our specific choice 
of GMSB models. It is indeed a remarkable coincidence 
that theoretically-motivated supersymmetric and grav- 
itino mass scales, the physical sizes of collider detectors, 
and gravitino cosmology all conspire to allow this corre- 
spondence. 

Note also that we find that the number of grav- 
itino events produced during the initial run of the LHC 
(center-of-mass energy of 7 TeV and integrated luminos- 
ity of 1 fb~^) may be comparable to that produced during 
an extended run of the Tevatron (center-of-mass energy 
of 2 TeV and integrated luminosity of 20 fb^^). This is 
true in regions of parameter space where large numbers of 
signals are expected. However, the higher center-of-mass 
energy of the LHC allows it to access regions of parame- 
ter space where rnsusY is larger; this is especially evident 
in the neutralino-NLSP scenario, as can be seen by com- 
paring the left and middle panels in both Figs. [2]and|3l 

If the distribution of mother-particle decay lengths can 
be measured with sufficient accuracy along with the total 
signal rate, then it may be possible to gain some infor- 
mation on the masses of the mother particles and the 
gravitino. To do so, it will be important to understand 
the distribution of energies and speeds with which mother 
particles are produced, since this will directly affect the 
distribution of decay lengths via dilation of the mother- 
particle lifetimes. In Fig. |6l we show some examples of 
probability distribution functions for the speed /? and the 
Lorentz factor 7 of mother particles that decay to grav- 
itinos, for various collider scenarios. 



V. CONCLUSIONS 

Light gravitinos in the mass range eV to MeV appear in 
GMSB models that naturally avoid flavor violation. We 
have examined the decay of supersymmetric particles to 
light gravitinos at colliders such as the Tevatron and the 
LHC. These decays will give rise to dramatic signatures, 
such as prompt di-photons or non-prompt photons, if the 
NLSP is a neutralino, or kinked charged tracks or heavy 
metastable charged particles, if the NLSP is a stau (or 
some other charged particle). We find large regions of the 
gravitino-mass-NLSP-mass parameter space in which the 
rate for such events may be appreciable at the Tevatron 
and LHC and which are consistent with current null su- 
persymmetry searches. 

Given that m^, ^ tonlsp for these events, the decay 
kinematics of individual events cannot be used to deter- 
mine the gravitino mass. However, the event rate and the 
distribution of decay locations may be used to narrow the 
range of NLSP and gravitino masses. Information about 
the nature of the NLSP may also be gleaned from the 
Standard Model decay products. 

One of the attractions of supersymmetry has been 
its ability to provide a natural candidate for the cold 
dark matter required by a wealth of cosmological ob- 
servations. Unfortunately, despite being well-motivated 
in GMSB models, the canonical light-gravitino scenario 
does not provide a natural cold-dark-matter candi- 
date. Nevertheless, this canonical scenario does allow 
gravitinos with masses ttiq < 30 eV that compose a 
fraction of the total dark matter, as determined by 
current astrophysical constraints on the relic abun- 
dance and small-scale structure. Given that upcoming 
structure formation observations are expected to probe 



9 



Tevatron (pp) - 2 TeV - 20 fb"' 




10'' 10" 10' 

mg (keV) 



650 
600 
550 
500 
450 
400 
350 



LHC (pp) - 7 TeV - 1 fb"'' 



to 

A 




1 


o 


X j 














/io - 










i 


■ 


i 


A < 80 TeV 



10'' 10" 10' 

mg (keV) 



650 
600 
550 
500 
450 
400 
350 



LHC (pp)-14TeV-10fb"'' 









'1 




to 


/r 


10 


7' 


e 


A 1 
X / 




\' / / 

/ / 




















If " 


- 10= ■—- 








JL-;'--- 10* 










. - - 10" ^ 


i 

^ 








A < 80 TeV 



prompt di-photon 
events 



SUSY events 
^ ^decav ^ 



10"' 10" 10' 

mg{keV) 



FIG. 2. Contour plots over the — rriNLSP parameter space showing the expected number of prompt di-photon events (soUd 
red) in a model with a neutralino NLSP, for the three collider scenarios (indicated at the top of each plot) of interest. The 
total number of supersymmetric events (dash-dotted black) and the average decay length ddccay (dashed black) expected at 
each point in the parameter space are also indicated by contours. The region in the upper-left corner is disallowed if we require 
A < A/mcss, while the region at the bottom is ruled out by LEP (using the conservative constraints discussed in the text). Note 
that hundreds of signal events may occur at the Tevatron with 20 fb~^ and at the 7 TeV LHC with 1 fb~^ if mg, < tens of eV. 
Observation of such a number of events would suggest that the canonical thermal-production scenario is correct, and that light 
gravitinos compose only a fraction of the dark matter. Also, note that a larger fraction_of neutralinos instead decay to ZG as 
the neutralino mass increases (and that there is some small fraction of decays to e'^e~G and hG). 
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FIG. 3. The same as Fig. (2] but for non-prompt photon events. Note that hundreds of signal events may occur at the Tevatron 
with 20 fb~^ and at the 7 TeV LHC with 1 fb~^ if tens of eV < mg < few keV. Observation of this number of events would 
suggest that a non-standard cosmology and gravitino thermal history cooled relic gravitinos, and, for the top part of this mass 
range, also diluted the relic density. 



hot-dark-matter masses as low as ttiq ~ eV, detection 
of a gravitino in the mass range eV < uIq < 30 eV via 
prompt signals at colliders vifould have implications 
for future small-scale-structure measurements. And 
although masses ttIq > 30 eV are disfavored, they may 
still be possible if the pre-BBN history was different 
than in the canonical scenario. Detection of gravitinos 
in this mass range via non-prompt and metastable 
signals at colliders would thus have serious implications 
for early-Universe cosmology, and may provide some 
insight into the reheating and inflationary eras. And 
who knows? There may indeed be new early-Universe 
physics that results in a gravitino that has the right 
cosmological abundance and temperature to be the dark 
matter. 
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FIG. 4. Contour plots showing the expected number of non-prompt lepton events in a model with a stau NLSP. Note that hun- 
dreds of signal events may occur at the Tevatron with 20 fb~^ and at the 7 TeV LHC with 1 fb~^ if tens of eV < tUq < few keV. 
Observation of this number of events would suggest that a non-standard cosmology and gravitino thermal history cooled relic 
gravitinos, and, for the top part of this mass range, also diluted the relic density. 
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